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ABSTRACT 

In the present paper we implement a non-trivial and renormalizable extension of 
the NJL model. We discuss the advantages and shortcomings of this extended model 
compared to a usual effective Pauli-Villars regularized version. We show that both 
versions become equivalent in the case of a large cutoff. Various relevant mesonic 
observables are calculated and compared. 



1. Introduction 



The Nambu and Jona-Lasinio (NJL) model Jljand its extensions have received much attention 
in low and medium energy hadronic physics [|2|,|3|. Because of its four fermion interaction it is 
non-renormalizable in the weak-coupling expansion for d > 2 dimensions, the reason for the model 
usually being treated with a cutoff A introduced to regularize the appearing ultraviolet (UV) 
divergencies. Although the original theory is non-renormalizable in perturbation theory, it becomes 
renormalizable in the mean-field expansion also for d > 2 |,|,|. However, in contrast to d < 4 
where the NJL model represents a perfect renormalizable field theory it is supposed to collapse 
for d = 4 to a trivial theory of non- interacting bosons [^|— p~0|] . Therefore in order to prevent the 
collapse, a cutoff A has also to be retained in the " renormalized" theory. The scale of this cutoff 
may in principle be deduced from an underlying non-trivial theory, in our case presumably QCD. 
If a finite A could be chosen large enough to extend the integrations of all finite integrals to infinity, 
then the cutoff appears only implicitly within the couplings which are prevented from being driven 
to zero. Formally this may be achieved by augmenting the model by bosonic kinetic terms and 
quartic self-couplings, capable to absorb the cutoff dependence of the coupling constants. This 
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procedure results in a renormalizable and non-trivial field theory for d = 4 dimensions which 
corresponds to a linear sigma model with quarks pj| , but with fixed bosonic self-couplings. In 
order to distinguish this model from the familiar cutoff or regularized NJL we will call it simply 
renormalized version in the following. 

This non-trivial extension of the NJL is motivated by the observation that some observables 
related to finite integrals require an infinite or at least a very large cutoff, most prominent example 
being the anomalous pion decay 7r° — > 77 [O . Due to the underlying symmetries of the model, it 
is clear that only processes which are sensitive to large momenta are sizeably influenced. Therefore 
many low-energy quantities, especially of course those used to fix the parameters of the model as 
e.g. the pion mass, are essentially the same as in the regularized version. 

The renormalized version has all the positive features known to renormalizable theories, but 
suffers from the occurrence of Landau ghosts, a well known problem related to Lagrangians without 
asymptotic freedom 

A recent work |l3| presents a renormalizable extension of the NJL-model by including the quark 
interaction generated by one gluon exchange, which simultaneously screens out the unphysical 
ghosts. The quark self-energy becomes momentum dependent with the appropriate asymptotic 
behavior, in contrast to the constant value obtained in the original NJL model. We consider our 
approach to be the "minimal" non-trivial renormalization program which still keeps the simple 
local structure of the original NJL Lagrangian and we think it worthwhile to analyse the results 
obtained from this approach, due to its simplicity. 

The aim of this paper is twofold. In the first part, which is of formal character, we show how to 
implement the original NJL Lagrangian to render it a non trivial renormalizable theory. Secondly 
we apply this Lagrangian to calculate relevant observables of the SU(2) flavor case and make a 
comparative study with the Pauli-Villars regularized version. We intend in this way to get a better 
understanding of the advantages and shortcomings of the two possible descriptions of the model 
in leading 1/Nc order. 



2. Mean field expansion of the NJL model and triviality 

In this section we briefly recapitulate the main features of the renormalization procedure for 
the NJL model using the mean-field expansion. Following Eguchi (J] we isolate the UV singulari- 
ties, but we consider also finite contributions which enter the renormalization scheme in order to 
demonstrate the equivalence to the procedure presented by Guralnik and Tamvakis S . In addition 
we allow the symmetry to be broken explicitly by a current quark mass. Finally we discuss the 
issue of triviality and the introduction of a cutoff A preventing the model from the collapse. The 
following derivation applies to SU(2) with N@ colors, an extension to SU(3) is straightforward. 

Starting point is the NJL Lagrangian with a local four-quark interaction 

L = fti-fd^ - m Q )q + ^ [(qq) 2 + {qi l5 rq) 2 } , (I) 

where we consider the isospin symmetric limit rh u = rhd = m. To allow for wave function renor- 
malizations later on and also in order to trace the Nc orders, it is convenient to replace Go = So/Mo 
with g 2 ~ 1/Nc- The subscript zero denotes bare (infinite) quantities everywhere. With boson 
fields introduced in the standard way the Lagrangian becomes 

L = q[i^d lJj - m - 50(00 + i^t-k )] q - y (00 + ^0) 

L = q[i^d„ - g (<7 + 275x71-0)] q - tt(°o + ^o) + — ™oVo , (2) 

* So 
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where the latter representation is obtained by shifting the scalar field <r — > °o — mo /go and a 
term independent of the dynamical fields is omitted. Integrations over the fermion fields q and q 
may now be performed in the path integral such that the resulting effective Lagrangian collects 
the corresponding trace log contribution 

2 2 

L = -iTrln - g (a + ij 5 Tir )} - ^(<Tq + 7T 2 ,) + ^W . (3) 

1 go 

Expecting the scalar field to possess a nonvanishing vacuum expectation value we expand cr = 
m/g + ct'q, where m represents the (finite) constituent mass 



oo 2 2 

L = i V -Tr - m)- 1 M + i 75 ttv )} " - ^K* + tt 2 ) - ^(m - m )a' . (4) 

„-i n 1 9o 



To evaluate this sum is now quite straightforward. The terms for n = 1, ... ,4 contain UV diver- 
gencies showing up as 



A L_ . = . f d A q 

(2^) 4 g 2 -m 2 ' l J (2^) 4 (q 2 - m 2 ) 2 



quadratically and logarithmically divergent integrals respectively. The effective Lagrangian is then 
given by 



\{-lNcglli og ){d^a' + ^tto^tto) - 

Z 1Z7T 



~(Mo - 8^c.9oVa<i)(^o 2 + *o) - ^m 2 (~4N c g 2 I log y 2 + ... 



+8N c goIio a m<7 Q (o' 2 + Tr 2 ,) + 2N c gUio a {<y'o + "l? + ■■■ 

o' Q . (6) 



n 2 

-^(to - m ) - 8N c gomI quad 

go 



The dots after the first two lines denote finite higher derivative terms of g 2 o 2 and g 2 ^ proportional 
to not explicitely shown, but a finite kinetic term for the scalars is kept and finally leads to 
different wave-function renormalizations a' = Z a 1 ^ 2 <t'q and tv = Z n 1 ^ 2 ttq. The dots in the third 
line indicate finite higher order terms also not shown. The rcnormalizcd parameters may now be 
introduced as follows 

4N c9 2 I log (7) 
iNcglhog (8) 
Z'gi(Iquad + 2m 2 I log ) (9) 

Cgllquad- (10) 

Together with the gap-equation 

-2 (to - mo) = SN c mI quad , (11) 
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which makes the linear term in a' vanish, we obtain the renormalized Lagrangian in its final form 



(12) 



The current quark mass has disappeared and it is noticed that the remaining parameters are related 
according to eqs. (7-10) 
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(13) 



such that the renormalized model is characterized by three parameters (m, g-n,n-n) as the regularized 
one (G, m, m). From these equations using the gap-equation (|ll]) we may also define a renormalized 
four-fermion coupling G 



rn 



9l 



9o 



m 
Go 



m 
G 



(14) 



by reintroduction of the renormalized (physical) current quark mass. Although this relation is 
beyond the scope of the renormalized model it will nevertheless prove useful for the evaluation of 
the quark condensate. 

The Nc orders are now carried by the couplings g 2 . ~ g 2 ~ 1/Nc- For practical calculations we 
consider only the leading order Nc for each process. From the quadratic terms in the Lagrangian 
we read off the renormalized meson propagators of order N c 



a;V)=p 2 -/4 



*N c g 2 a 



(p 2 - 4m 2 )Z (p 2 ) 



V 



24tt 2 



A" V) =V-A- ±N c gip 2 Z (p 2 ) , 



(15) 



where the momentum dependent terms due to the finite higher derivative terms not explicitely 
shown in eq.([l2|) are contained in the finite function 



Z ( P 2 ) 



1 



16tt 2 
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8^ 



dz£n 



1 



P 



z(l-z) 



1 + 



4m 2 



FT 



arsinh 



4 m - 

V 1 



1 - 



1 arcsin 
K cosh 



|P 2 | 

4m 2 
Ant 2 



Ini- 



p 2 < 
< p 2 < 4m 2 
4m 2 < p 2 



(16) 



which possesses different branches. The physical meson masses are defined as usual via the poles 
of the propagators 



ml = Am 2 +nl+ 4,N c g 2 Jml - 4m 2 )Z (m 2 ) 
ml =/i 2 + 47V c . 9 2 m 2 Z (m 2 ). 

In the chiral limit (/i 2 = 0) we obtain m^ — and m a = 2m. 



(17) 
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Similarly, also in the chiral limit the leading Nq 3- and 4-boson vertex functions at zero momenta 
are obtained from the cubic and quartic self-couplings in accordance with Guralnik and Tamvakis 
H, who derive these results from the corresponding Ward identities. 

In the following we want to discuss the triviality of the NJL model which is suggested by lattice 



calculations |10|. In the mean-field expansion it follows immediately from (8) 

gl = -(4N c I log )-\ (18) 

namely in the continuum limit the couplings and hence also g a are driven to zero rendering 
the Lagrangian @ a theory of non-interacting mesons. This is caused by the mesonic kinetic 
terms and quartic self couplings in ([l2]) being created purely by radiative corrections: they were 
not present in the original Lagrangian. To avoid the collapse of the model, g^ must be kept 
fixed at some finite value. For that purpose a cutoff A has to be introduced in order to keep the 
logarithmically divergent integral in ( p"8| ) finite (note that the quadratic divergence has already 
disappeared in the renormalized parameters). In the continuum limit g 2 ~ (AiT/Nc)/£n(k/m) 
tends to zero logarithmically, and a finite coupling g n requires also a finite A, in fact, in order to 
reproduce a reasonable coupling strength a rather low cutoff of the order of lGeV is needed in 
contrast to the situation in QED where the collaps is prevented by a cutoff located way above all 
physical energies of interest. Nevertheless, if it were possible to choose A large enough such that 
all finite integrals may be evaluated in the continuum limit then the cutoff would disappear from 
the theory being only implicitely contained in the coupling g n which is kept finite |L4]]. Exactly 
this is achieved by adding mesonic kinetic terms and quartic self-couplings to the model which are 
capable to absorb the troublesome radiative terms in ( |l2| ) leading to a non-trivial renormalizable 
extension of the NJL discussed in the following section. Mesonic properties calculated in the two 
versions of the model are then presented in section 4. 



3. Non-trivial extension of the NJL model 



We have seen in the previous section that the triviality of the NJL model is connected with the 
fact that the mesonic kinetic and interaction terms are created purely by radiative corrections. 
In fact triviality may be avoided by adding these contributions to the Lagrangian (||) from the 
beginning 



2 

ii 2 \ ,, 2 

PO/ 2 2\ A / 2 , 2\2 i Mo 



-iTrin [i'fd^ - g {a + «75™ )] + ^(S^o + d^nd^o) 



* 1 9o 

Of course this may lead beyond the NJL model, we will comment on this later. Repeating the 
steps which lead to eq. (^) we find the renormalized parameters as 

Z- 1 = § = fl ANcglh 09 - ^ (20) 

2 

Z~ X = f| = /O - ±Nc9tllo g (21) 



U 6A 

^ = Mo - %N c gl{I q uad + 2m 2 I log ) + -^-m 2 (22) 



u 2A 

^ = Mo - 8N c g 2 I quad + -^m 2 (23) 

Z 7T 9q 
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9Z % 
together with the gap-equation 



"2 = -J — INcIlog (24) 



2 ^ (m - m ) - 8N c mI quad + -^to j = . (25) 
5o 9o 



The renormalized Lagrangian in the shifted scalar fields becomes 



.^V'Q^' 2 + gW) ^(g^ + gWf +... (26) 



formally identical to the bosonized NJL Lagrangian (|l2j) if we choose A = 1. In general the model 
parameters are now related 

_1_ = 1_ _ Nc_ /4 = 4Am 2 + A 

gl gl 6tt 2 ' gl g\ 

quite similar to (|l^) and also the relation ( |l4| ) remains unchanged. We want to emphasize however, 
that this Lagrangian does not represent the trivial NJL model and constitutes instead a different 
non-trivial theory. As a consequence the couplings g v and g a do no longer vanish in the continuum 
limit. The numerical results of this model with A = 1 are of course identical to those obtained from 
the NJL model ( ^2|) with g^ and g a kept fixed at some finite values as discussed in the preceeding 
section. 

Concluding, there seems to be two options to treat the problem of triviality which appears in 
the NJL model: 

(i) A cutoff A is retained to prevent the model from the collapse. Because numerically the cutoff 
is of the order of 1 GeV only, it has to be kept also in all finite integrals. 

(ii) The NJL model is augmented by kinetic terms and mesonic self-interactions. This results in 
the linear sigma model coupled to quarks and constitutes a perfect non-trivial field theory. 

The latter theory contains one additional parameter A. For A = I (ii) gives the same results as the 
conventional NJL in case of a large cutoff. In principle the assumption A = 1 may be tested in mr 
scattering, of course not in the leading chiral order which is fixed by a famous low energy theorem 
|l6| ], but in the next to leading orders (subsection 4.2). Many other meson properties are quite 
independent of this parameter. In the following section we compare some mesonic observables 
calculated in the two versions of the NJL model. 



4. Meson properties 

In this section the parameters of both versions of the model are fixed. We show that in the 
renormalized version the chiral expansion becomes quite simple and we discuss the issue of the 
additional parameter A appearing in the linear sigma model with quarks. Finally we are going to 
calculate several mesonic observables. 



G 



1. Determination of the model parameters 



For both versions of the model we used two sets of model parameters : one with a low constituent 
quark mass fixed at to — 210 MeV, and the other with the constituent quark mass fixed at to = 350 
MeV, values which are used widely in the literature. The remaining model parameters are adjusted 
to reproduce the pion decay constant / w = 93.3 MeV and the pion mass m v — 139 MeV. As a 
result, one obtains for the Pauli-Villars regularized version a large dimensionless ratio A/to ~ 6 
in the small constituent quark mass case, and a small one for the large mass case, A/to ~ 2, see 
Table I. As will be discussed in section 4.2, chiral expansions of the renormalized and regularized 
versions coincide in the A/to — » oo limit, showing that this ratio is a measure for the deviations in 
the two models. On the other hand, in the renormalized version the pion decay constant is 



9-irqq o 2 

g£ mi 



(28) 



and the pion mass is given by (17) 
quark couplings 



Here and in the following we use the pion quark and sigma 



g-* q = g- 2 - 4N C [Z (m*) + m£Z$(m*)] 

9^ q = 9^ - 4A C [Z (mD + (m» - 4to 2 )^(to^)] 



(29) 



for abbreviation. Thus, f v and m v fix the parameters g^ and fi^ listed also in Table I. Furthermore, 
for the evaluation of the quark condensate according to ( |l4| ) 



TO P" 

to < qq >= — (to — m)— = — m(m — to)— ^ 

G gi 



(30) 



a value for the current quark mass has to be adopted which strictly speaking is not a parameter 
of the renormalized version of the model. The standard value of to = 7.5 MeV leads to the result 
quoted in Table II. 



TABLE I. Parameters of the regularized and renormalized versions of the NJL model. 





regularized 


model 




renormalized 


model 




model 


m = 350 


(210) 


MeV 


m = 350 


(210) 


MeV 


parameters 


G = 17.6 


(5.11) 


GeV" 2 


g„ = 3.752 


(2.250) 






m = 8.5 


(4.1) 


MeV 


H-k = 141 


(141) 


MeV 


related 


A = 769 


(1190) 


MeV 


g a = 7.006 


(2.610) 




parameters 








/i CT = 1333 


(513.8) 


MeV 
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2. Chiral expansion and tttj scattering 



In the renormalized version of the NJL model the chiral expansion for constituent quark mass, 
pion decay constant, pion mass and quark condensate become quite simple: 



rn 



J 7T 



4m 



2 
At 



m * N cm T] 



2m 



mZ = m„ 



m < qq >- 



o 2 

Ncm„ 



4m 12 ^ 2/t 



J 7T 71 



rn 

o 

rn 



4m 



(31) 



(32) 



(33) 



(34) 



These formulas agree with those obtained in the Pauli-Villars regularized model [17]] for large 
cutoff A 3> m. It is noticed that the current quark mass appears only in the expression for the 
condensate. 

Similarly, the 7T7t scattering amplitude (box-diagram + sigma exchange) is obtained as 



A(s,t,u) = 



o2 
At 



1 - 



N c m 



N c m 
^2" i ^4 



o 2 Ncm 
s-m„ 7T 



2m, 



4m / ff 



4m /„ 



2m, 



Nr 



24tt2/, 



s(u + t) — ut — 2m w 



(35) 



As mentioned already, the leading term is fixed by a low energy theorem |l6[ ] and is therefore 
independent of the strength A of the quartic self-interaction. However, in the next to leading order 
there appears a term in addition to those obtained in the regularized model with large cutoff, 
which is effective for A ^ 1. 

In refs. Jlj],|| it was shown that compatibility with chiral perturbation theory (ChPT) requires 
a small constituent quark mass of the order of m ~ 250 MeV. For such a small quark mass the 
cutoff is of minor importance and the terms in ( |35| ) which survive for A = 1 fit the ChPT nn 
scattering threshold parameters fl8|| in the renormalized version as well. We may conclude that 
the comparison with ChPT requires a value of A close to 1. Of course, this conclusion is correct 
only if we accept a small constituent quark mass. 



3. Sigma meson properties 

The sigma meson mass can be evaluated using eq. (|T^) . In order to obtain this mass, we made 
use of the real part of the propagator only, in both, the regularized and the renormalized version. 
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The imaginary part related to the decays into qq pairs is very small and can be neglected [ |19| . The 
decay width of a — ► ttt: is : 



3a/ m 2 — Am 2 2 



87rm2 



r 

J (J7T7T 



(36) 



All the expressions for the Pauli-Villars regularized NJL model used here and in the following sec- 
tions may be found in po| and are not repeated here. The amplitude f avv reads in the renormalized 
version 



/cttttt = l6mN c ga-q q gl qq 



Zo(ml)+ 2m2 \ ml I 3 (ml,ml,ml) 



h{p 2 ,p\,p\)=i 



d 4 q 



(2ir) 4 (q 2 — m 2 )[{q — p\) 2 — m 2 ] [{q — pi) 2 — m? 



(37) 



with p 2 = (pi -P2) 2 . 

As is noticed from Table II, the scalar decay width into two pions depends most sensitively on 
the two versions of the model, specially of course in the large mass (low cutoff) case. This difference 
persists also in the chiral limit. 



4- Pion charge formfactor 



In the renormalized version, the pion electromagnetic formfactor is given by 



F v (p 2 ) = 2N c gl Q 



1 



AN c gl 



+ ml (Z' (ml) - h{p 2 , ml, ml)) 



(38) 



and the corresponding pion electromagnetic radius is defined as usual. In Fig. 1 the electromagnetic 
pion formfactor is shown in the space-like region only, since in the time-like region the role of vector 
mesons, which are not considered here, is very important. For the same reason the charge radius 
turns out too small compared to its experimental value, see Table II. We see however that the 
renormalizable version yields a value close to the chiral limit result 



3A^ 
^ 2 f 2 



(0.59/m) 2 



(39) 



whereas the regularized version tends to decrease this value further with decreasing ratio A/m. 
From Fig. 1 we see that the formfactors in the two versions of the model are very similar for the 
lower mass case, but for the larger mass case the renormalized version improves on the results, 
whereas the regularized version does not yield a satisfactory fit as noticed already in 02] . 
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FIG. 1. Pion formfactor in the space- like region. The results of the renormalized (solid) and regularized 
(dashed) models for small and large constituent quark masses are compared to experimental data 

5. Anomalous n — * 77 decay 

The formfactor associated with the anomalous process no — ► 7*7, with one of the photons being 
off shell is 



8N c e- 



9irq q ml 3 (0,p ,771,,) 



(40) 
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With both photons on-shell one obtains the anomalous pion decay 7To 
renormalized version, given by the following expression: 



+77 



//to O 



N c e 2 g„ 



: arctan ■ 



77 analytically in the 



(41) 



In the chiral limit, this reduces to the well-known result [E2l 



r„ 



7T U 



+77 



N c e 2 
64tt V12tt 2 / ot 



(42) 



The model independent amplitude for the anomalous ir° — > 77 decay in the chiral limit is ob- 
tained exactly only in the renormalized version. The regularized version renders the decay width 
strongly cutoff dependent (see Table II). This formal result is one of the obvious advantages of the 
renormalized model. In this context we should mention that in the spirit of regularized models a 
consistent treatment of anomalous processes has been forwarded in [23)| , where the quark-loop is 
regulated dynamically by the intrinsic non-locality of the quark-meson interaction. 

The 7T° — > 7*7 transition formfactor as it is obtained in the renormalized version is plotted in 
Fig. 2. We do not compare to the results of the regularized version, because it fails already at 
the photon point (Table II). Obviously here the lower mass case leads to a much better fit to 
experiment. The slope of the formfactor at the origin is obtained as a„o — 0.035to^ o which has to 
be compared to the empirical value a^o = 0.0326 ± 0.0026to^ q , p4j . 



TABLE II. Some meson properties calculated in the regularized and renormalized versions of the NJL 
model are compared to experimental data [ piU^ ]. The asterisks indicate quantities which served as input 
to determine the model parameters. Results are calculated for a constituent quark mass of m = 350 MeV 
and m = 210 MeV (in brakets). 





regularized 


model 


renormalized 


model 


experiment 


u 


[MeV] 


93.3* 




93.3* 




93.3 




[MeV] 


139* 




139* 




139 




[MeV] 


705 


(434) 


705 


(433) 


~ 700 


r 


[MeV] 


647 


(336) 


457 


(307) 


~ 700 


L 7T° — 


[eV] 


3.7 


(7.2) 


7.6 


(8.3) 


7.7 ±0.6 




[fin] 


0.49 


(0.59) 


0.58 


(0.59) 


0.678 ±0.012 


<qq> 


[MeV 3 ] 


-265 3 


(-340 3 ) 


-283 3 


(-281 3 ) 


-(270 ±20) 3 
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FIG. 2. Transition formfactor for the anomalous process ty° — > 7*7 in the space-like region as calculated 
in the renormalized version with different constituent quark masses. The data are from Ref. |24j . 

For the larger A/to ratio the results in Table II are quite similar in the two versions as expected. 
For the larger constituent mass case, however, the results for the anomalous pion decay and the 
pion electromagnetic radius are substantially improved in the renormalized version. 



5. Conclusions 



We have considered a non-trivial and renormalizable extension of the NJL model obtained by 
adding the necessary counter terms, namely the mesonic kinetic energies and quartic mesonic 
interactions to the originial Lagrangian. This amounts to a linear sigma model coupled to quarks 
with fixed strength of the quartic mesonic self-interactions. It was shown that this version of the 
model coincides with its familiar regularized versions provided the cutoff is large enough. 

We presented a comparative and quantitative analysis of the most relevant mesonic observables 
with an effective Pauli-Villars regularized version in leading 1/Nq order. We found that in the 
renormalized version of the model the electromagnetic pion decay is always (independently of other 
parameters) in agreement with experiment in contrast to results obtained in regularized versions 
(usually > 30% off). For the first time a reasonable description of the corresponding transition 
formfactor is obtained. In the regularized versions the quality of the results depends crucially 
on the ratio A/m where A is the four-dimensional cutoff and to is the constituent quark mass. 
The larger the ratio, the better the agreement with phenomenology. This is in a way reflected in 
the renormalized version, where the agreement is systematically better when compared with the 
regularized version for the same m value. We noticed that the quantitative results become quite 
similar in the two versions for ratios A/m > 6. 

Further applications of the extended non-trivial and renormalizable NJL model are the study 
of the SU(3) flavor case and the inclusion of vectors and axialvectors mesons. One of the most 
interesting issues related to the SU(3) sector is connected with the kaon decay constant which 
persistently turns out to be too small (fx — fir) unless one is willing to accept a small nonstrange 
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constituent quark mass of the order m ~ 210MeV which brings along its own difficulties (small 
scalar meson masses, low qq threshold, etc.). Unfortunately the renormalizable version presented 
here is not capable to resolve this problem, the solution has to be sought elsewhere. 

We are grateful to O. A. Battistel for critical remarks on regularization prescriptions p7fl . 
H.W. would like to thank G. Holzwarth for clarifying comments on the renormalization of 
the model. The present research has been supported partly by JNICT, Portugal (Contract 
PRAXIS/4/4.1/BCC/2753 and PCERN/S/FIS/1162/97) and by CNPq-Brazil. 
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